Carbon clusters near the crossover to fullerene stability 
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The thermodynamic stabihty of structural isomers of C24, C26, C28 and C32, including fuUerenes, is studied using 
density functional and quantum Monte Carlo methods. The energetic ordering of the different isomers depends 
sensitively on the treatment of electron correlation. Fixed- node diffusion quantum Monte Carlo calculations 
predict that a C24 isomer is the smallest stable graphitic fragment and that the smallest stable fuUerenes are the 
C26 and C28 clusters with C2v and Td symmetry, respectively. These results support proposals that a C28 solid 
could be synthesized by cluster deposition. 
PACS: 61.46. +w, 36.40.-C, 71.45. Nt 



Since the discovery of the fullerene Cgo , the study 
of carbon clusters has revealed a rich variety of physical 
and chemical properties. Fullerene clusters may now be 
synthesised in macroscopic quantities, but despite many 
experimental and theoretical advances the detailed en- 
ergetics of these systems are not yet fully understood. 
The question "which is the smallest stable fullerene?" 
remains both interesting and contentious due to the sen- 
sitivity of cluster formation to experimental conditions 
and the challenges posed to theoretical methods by sys- 
tem size and the high accuracy required. In this Letter 
we report very accurate calculations of the relative ener- 
gies of C24, C26, C28 and C32 clusters, and identify the 
smallest stable fuUerenes. 

The number of low-energy candidate structures can be 
large, even for quite small clusters, precluding exhaus- 
tive theoretical searches with highly accurate but com- 
putationally expensive methods. In practice, a hierar- 
chy of methods of increasing accuracy and computational 
cost must be used. The first step is to select candidate 
structural isomers via empirical methods based on bond 
counting and geometric "rules" such as "minimize the 
number of adjacent pentagons" j2j. Quantum mechani- 
cal calculations based on tight-binding and density func- 
tional theory (DFT) methods can then used to refine the 
selection. To finally establish the energetic ordering of 
different isomers, highly accurate calculations must be 
performed. Quantum chemical methods, such as coupled 
cluster (CC) calculations ||], are potentially highly accu- 
rate, but are severely limited by the size of basis set that 
is computationally affordable in these systems. Quantum 
Monte Carlo (QMC) methods give an accurate treatment 
of electron correlation which, combined with an absence 
of basis set error, favorable scaling with system size and 
suitability for parallel computation, renders them ideal 
for these studies. QMC calculations have reproduced 
experimental binding energies of small hydrocarbons to 
within 1% Using the techniques described below we 
have calculated the cohesive energy of bulk diamond, ob- 
taining values of 7.36(1) and 7.46(1) eV per atom in vari- 
ational Monte Carlo (VMC) and diffusion Monte Carlo 



(DMC), respectively, which are in very good agreement 
with the experimental value of 7.37 eV. 

Carbon clusters are particularly challenging to model 
accurately due to the wide range of geometries and the 
occurrence of single, double, and triple bonds. These dif- 
ferences result in a non-cancelation of errors in relative 
energies, exaggerating any errors due to approximations 
involved in electronic structure methods. Despite these 
potential difficulties, carbon clusters have been exten- 
sively studied using methods such as tight-binding, den- 
sity functional, quantum chemical and QMC The 
need for high accuracy calculations with a sophisticated 
treatment of electron correlation has been clearly illus- 
trated by several previous studies. Grossman et al. ^ 
have performed diffusion Monte Carlo calculations for 
C20 clusters, finding that the fullerene is not energetically 
stable. A DFT study of C20 isomers showed that the local 
density approximation (LDA) [|j and the BLYP gradient- 
corrected functional gave different energy orderings 
for the sheet, bowl and fullerene structures, with nei- 
ther agreeing with that of DMC Q. Jensen et al. |l|] 
made calculations for the monocyclic ring and fullerene 
isomers of C24 which demonstrated significant differences 
between the predictions of the LDA, gradient corrected 
and hybrid density functionals. These authors also per- 
formed second-order M0ller-Plesset and CC calculations, 
but with a limited basis set, concluding that the fullerene 
is lower in energy than the ring. Raghavachari et al. |l^ ] 
and Martin et al. |l^ studied seven isomers of C24 using 
DFT, but obtained conflicting energetic orderings. 

For clusters containing between 20 and 32 atoms, three 
classes of isomer are energetically competitive: fuUerenes, 
planar or near-planar sheets and bowls, and monocyclic 
rings. The smallest possible fullerene, defined as a closed 
cage containing only pentagonal and hexagonal faces , 
consists of 20 atoms. However, the smallest fuUerenes 
most commonly identified by time of flight and mass 
spectroscopy measurements are the C30 and C32 clus- 
ters ]l5|-p7|. Rings are found to dominate up to ap- 
proximately 28 carbon atoms under typical experimental 
conditions, and fuUerenes are mostly observed for larger 
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clusters, although other structures are also present (see 
for example, Ref. In this work we present a QMC 

study of five isomers of C24, three of C26 and C28, and 
two of C32, thereby covering the range of masses from 
where the fuUerene is clearly predicted to be unstable to 
where a fuUerene is clearly observed. This enables us to 
predict the smallest energetically stable fuUerene. 

We apply the diffusion quantum Monte Carlo method 
1^,^ in which the imaginary time Schrodinger equation 
is used to evolve an ensemble of electronic configurations 
towards the ground state. The "fixed node approxima- 
tion" is central to this method; the nodal surface of the 
exact fermionic wave function is approximated by that of 
a guiding wave function. Core electrons were modeled by 
an accurate norm-conserving pseudopotential and 
the non-local energy was evaluated stochastically within 
the locality approximation We used Slater- Jastrow 
guiding wave functions consisting of the product of a sum 
of Slater determinants of single-particle orbitals obtained 
from CRYSTAL95 p|l or Gaussian94 [|3) with a Jastrow 
correlation factor ]24| ]. Optimized uncontracted valence 
Gaussian basis sets of four s, four p and one d function 
were used to represent the single-particle orbitals. Jas- 
trow factors of up to 80 parameters were optimized using 
efficient variance minimization techniques |^,^, yield- 
ing 75-90% of the DMC correlation energy. 

We relaxed the structures by performing highly con- 
verged density functional calculations. The geometries 
were obtained from all-electron calculations ||2^ using the 
B3LYP hybrid functional ^ and Dunning's cc-pVDZ 
basis set |28| , which has been found to be an accurate and 
affordable combination . To assess the sensitiv- 

ity of the total energies to the geometries, we compared 
the energies of the fully relaxed ring and Dg fuUerene 
isomers of C24 (see Fig. |^) using the BLYP and B3LYP 
functionals. The functionals give significantly different 
energetic orderings, but the differences between the ge- 
ometries are small - less than 0.03 angstroms in bond 
lengths and 0.4 degrees in bond angles. The relative en- 
ergies of these structures changed by a maximum of 0.27 
eV for each of the functionals investigated. The relative 
energies are therefore rather insensitive to the functional 
used to obtain the geometries, but are more sensitive 
to the functional used to calculate the energies. These 
changes are small compared with the overall range of en- 
ergies, but some changes in the orderings of the isomers 
closest in energy could occur. 

We considered the following isomers of C24, as de- 
picted in Fig. |l|: a polyacetylenic monocyclic ring, a flat 
graphitic sheet, a bowl-shaped structure with one pen- 
tagon, a caged structure with a mixture of square, pen- 
tagonal and hexagonal faces, and a fuUerene. Other can- 
didate structures, such as bicyclic rings and a 3-pentagon 
bowl were excluded on the grounds that DFT calculations 
using several different functionals have shown them to be 
significantly higher in energy As well as DMC 



calculations we have also performed DFT calculations 
using the LDA, two gradient corrected functionals (PEE 
|l| and BLYP) and the B3LYP functional. The resuUs 
shown in Fig. |l| confirm that the treatment of electron 
correlation has a profound effect on the relative energies. 
All of the functionals give different energetic orderings, 
and none gives the same ordering as DMC. The graphitic 
sheet is placed lowest in energy by DMC, in agreement 
with each of the functionals except BLYP, which places 
the ring lowest in energy. The low energy of the C24 
graphitic sheet is expected because the structure accom- 
modates a large number (7) of hexagonal rings without 
significant strain. This structure is predicted to be the 
smallest stable graphitic fragment. Both DMC and the 
DFT approaches find the C24 fuUerene to be unstable. 

Three isomers of C26 were considered: a cumulenic 
monocyclic ring, a graphitic sheet with one pentagon and 
a fuUerene of C2V symmetry (Fig. ||). Few studies of the 
C26 fuUerene have been made, in part due to the high 
strains evident in its structure . Recently Torelli and 
Mitas have demonstrated the importance of using multi- 
determinant trial wave functions to describe aromaticity 
in 4N-I-2 carbon rings We have tested this for the 
C26 ring, using a 43 determinant trial wave function ob- 
tained from a CI singles-doubles calculation. The multi- 
determinant wave function gave a slightly lower DMC 
energy than the single determinant wave function, by 
approximately 0.5 eV, confirming that CI wave functions 
can have better nodal surfaces than HF wave functions. 
The ring and sheet-like isomers are close in energy, but 
the fuUerene is approximately 2.5 eV below these iso- 
mers and is therefore predicted to be the smallest stable 
fuUerene. Small changes in the geometries are highly un- 
likely to change this conclusion. 

Three C28 isomers were investigated: a monocyclic 
ring, a graphitic sheet and a fuUerene of Tj symmetry 
(Fig. Other bowl and sheet-like structures were ex- 
cluded on energetic grounds Spin-polarized DFT 
calculations show the ground state of the Td symmetry 
fuUerene to be a spin-polarized ^^2 state. DMC predicts 
that this spin-polarized fuUerene is the lowest energy iso- 
mer of C28, and this is supported by each of the func- 
tionals except BLYP. The spin-polarized fuUerene has 
four unpaired electrons and is therefore highly reactive. 
This property has been exploited in atom trapping exper- 
iments in which fuUerenes containing single four-valent 
atoms, C28M, have been prepared by laser vaporization 
of a graphite-M02 (M — Ti, Zr, Hf or U) composite rod 
p3[ . Our prediction that the fuUerene is the most stable 
isomer of C28 indicates that isolated fuUerenes might be 
readily produced. This would facilitate investigations of 
C28 fuUerene solids, which have been discussed but not 
yet produced |3^ , |3^ , although this route may be ham- 
pered by the chemical reactivity of the fuUerene. (A C36 
fuUerene solid has been reported |^.) 

Our DFT and DMC resuUs for C28 (Fig. ||) again high- 
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light a wide variation between different DFT functionals. 
The LDA and B3LYP functionals predict the same or- 
dering as DMC, but the PBE and BLYP functionals give 
different orderings. The DMC data strongly indicates 
that the Td fuUerene is the most stable C28 isomer at 
zero temperature. The fuUerene has the lowest DMC en- 
ergy in both spin-polarized and non spin-polarized calcu- 
lations, and is substantially more stable than the sheet 
and ring. Small changes in the geometries are therefore 
unlikely to change this ordering. 

Our DMC calculations for the C32 monocyclic ring and 
fuUerene show that the fuUerene is 8.4(4) eV per molecule 
lower in energy, which is consistent with the observation 
of a large abundance of C32 fuUerenes in a recent clus- 
ter experiment In Fig. ^ we plot the DMC binding 
energies per atom of all the ring and fuUerene structures 
considered. The binding energies of the fuUerenes rise 
much more rapidly with cluster size than those of the 
rings because of the large amount of strain in the smaller 
fuUerenes. The DMC binding energy of the C32 fuUerene 
is approximately 1 eV per atom less than the experimen- 
tal binding energy of Cgo- 

Our DFT and DMC results highlight several impor- 
tant trends in the relative performance of the different 
functionals. The overall quality of a functional for the 
clusters is best judged by the agreement with the DMC 
data for the overall shapes of the relative energy data of 
Figs. ^[^. The best agreement is given by the PBE and 
B3LYP functionals, with the LDA being slightly inferior 
and the BLYP functional being worst. The tendency of 
the LDA to favor structures of high average coordination 
and for the BLYP functional to favor structures of low 
average coordination is consistent with the results for C20 
reported by Grossman et al. Q . 

The final test of our predictions must lie with experi- 
ment. It is clear that the actual abundances of different 
clusters depend sensitively on experimental conditions. 
Analysis of the stability of clusters against fragmenta- 
tion, growth and other chemical reactions is complicated. 
One key issue is that the clusters are formed at tem- 
peratures of order lOOOK and therefore the vibrational 
contributions to the free energy can be significant. For- 
tunately, a simple picture emerges from computations of 
vibrational properties jl^ jl^,^ . FuUerenes are relatively 
rigid and have higher vibrational free energies than rings, 
which have many low-lying vibrational modes. Vibra- 
tional effects therefore tend to favor the ring isomers at 
high temperatures. However, according to our DMC cal- 
culations the C26 and C28 fuUerenes are several eV per 
cluster lower in energy than the other isomers, so that 
significant amounts of fuUerene could exist at the tem- 
peratures of formation. If thermodynamic stability alone 
were to determine which cluster sizes were observed then 
only the largest fuUerenes would ever be observed, but in 
a recent experiment the abundance of the C32 fuUerene 
was found to be greater than Cgo fisl] . There is more evi- 



dence that thermodynamic stability to rearrangements of 
clusters of a particular size are important in determining 
which isomers are observed. For example, in the experi- 
mental study of Ref . , fuUerenes were mostly observed 
for clusters containing more than about 30 carbon atoms, 
while for smaller clusters mostly rings were formed. This 
crossover is close to the critical size for fuUerene stability 
of 26-28 atoms predicted by our DMC calculations. 

In conclusion, performing accurate calculations of the 
relative energies of carbon clusters is a severe test of elec- 
tronic structure methods because of the widely differing 
geometries and the occurrence of single, double and triple 
bonds. In our DMC calculations for C24, the lowest en- 
ergy isomer is a graphitic sheet, which is expected to be 
the smallest stable graphitic fragment. We predict that 
the smallest energetically stable fuUerenes are the C2V 
symmetry C26 cluster and the spin polarized ^ A2 state 
of the Td symmetry C28 cluster. This prediction lends 
weight to recent proposals that a C28 solid could 
be synthesized by surface deposition of C28 fuUerenes. 

Financial support was provided by EPSRC (UK). Cal- 
culations were performed on the CRAY-T3E at the Uni- 
versity of Manchester and the Hitachi SR2201 located at 
the University of Cambridge HPCF. 



[1] H. W. Kroto, J. R. Heath, S. C. O'Brien, R. F. Curl, and 

R. E. Smalley, Nature (London) 318, 162 (1985). 
[2] H. W. Kroto, Nature (London) 329, 529 (1987). 
[3] J. Cizek, Adv. Chem. Phys. 14, 35 (1969). 
[4] J. C. Grossman, L. Mitas, and K. Raghavachari, Phys. 

Rev. Lett. 75, 3870 (1995); Erratum: Phys. Rev. Lett. 

76, 1006 (1995). 
[5] J. Cioslowski, Electronic Structure Calculations on 

FuUerenes and Their Derivatives (Oxford University 

Press, 1995). 
[6] G. E. Scuseria, Science 271, 942 (1996). 
[7] A. Van Orden and R. J. Saykally, Chem. Rev. 98, 2313 

(1998). 

[8] Y. Shiakhter et al, J. Chem. Phys. 110, 10725 (1999). 
[9] J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 
(1981). 

[10] A. D. Becke, Phys. Rev. A 38, 3098 (1988); C. Lee, 
W. Yang, and R. G. Parr, Phys. Rev. B 37, 785 (1988). 

[11] F. Jensen and H. Koch, J. Chem. Phys. 108, 3213 (1998). 

[12] K. Raghavachari et al, Chem. Phys. Lett. 220, 385 
(1994). 

[13] J. M. L. Martin, Chem. Phys. Lett. 255, 7 (1996). 

[14] P. W. Fowler and D. E. Manolopoulos, An Atlas of 

FuUerenes (Clarendon Press, Oxford, 1995). 
[15] H. Kietzmann et al, Phys. Rev. Lett. 81, 5378 (1998). 
[16] H. Handschuh et al, Phys. Rev. Lett. 74, 1095 (1995). 
[17] G. von Helden et al, Chem. Phys. Lett. 204, 15 (1993). 
[18] D. Ceperley, G. Chester, and M. Kalos, Phys. Rev. B 16, 



3 



3081, (1977). 

[19] B. L. Hammond, W. A. Lester, Jr., and P. ,J. Reynolds, 
Monte Carlo Methods in Ab Initio Quantum Chemistry, 
(World Scientific, Singapore, 1994). 

[20] N. TrouUier and Jose Lui's Martins, Phys. Rev. B 43, 
1993 (1991). 

[21] M. M. Hurley and P. A. Christiansen, J. Chem. Phys. 

86, 1069 (1987); B. L. Hammond, P. J. Reynolds, and 

W. A. Lester, Jr., J. Chem. Phys. 87, 1130 (1987); 

L. Mitas, E. L. Shirley, and D. M. Ceperley, J. Chem. 

Phys. 95, 3467 (1991). 
[22] R. Dovesi et al, CRYSTAL95 User's Manual (University 

of Torino, Torino, 1996). 
[23] M. J. Frisch et al, Gaussian 94 (Gaussian Inc., Pitts- 
burgh, PA, 1995). 
[24] A. J. WiUiamson et al., Phys. Rev. B 53, 9640 (1996). 
[25] C. J. Umrigar, K. G. Wilson, and J. W. Wilkins, Phys. 

Rev. Lett. 60, 1719, (1988). 
[26] P. R. C. Kent, R. J. Needs and G. Rajagopal, Phys. Rev. 

B 59, 12344 (1999). 
[27] A. D. Becke, J. Chem. Phys. 98, 5648 (1993). 
[28] T. H. Dunning Jr., J. Chem. Phys. 90, 1007 (1989). 
[29] J. M. L. Martin et al, Mol. Phys. 86, 1437 (1995). 
[30] J. M. L. Martin, Chem. Phys. Lett. 255, 1 (1996). 
[31] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 

Lett. 77, 3865 (1996). 
[32] T. Torelli and L. Mitas (unpublished). 
[33] T. Guo et al, Science 257, 1661 (1992). 
[34] J. Kim et al, J. Chem. Phys. 108, 2631 (1998). 
[35] C. Piscotti, T. Yarger, and A. Zettl, Nature (London) 

393, 771 (1998). 



> 

0) 



C 
<B 

> 
(D 




A BLYP 
▼ B3LYP 



Ring Sheet Bowl O,, Cage Dg Fullerene 

FIG. 1. The structures and energies of the C24 isomers 
given relative to the Dg fullerene. 
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The structures and energies of the C26 isomers. 
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FIG. 3. The structures and energies of the C28 isomers, 
given relative to the spin-polarized fullerene (see text). 
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FIG. 4. The DMC binding energies of the C24, C26, C2S, 
C32 ring and fullerene structures. The lines drawn arc for 
guidance only. Statistical error bars are smaller than the sym- 
bols. 
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